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! Abstract 

o 

We calculate the leading 1/iVy perturbative contributions to the 
polarized nonsinglet structure functions 53 and 55 to all orders in a s . 



The contributions from the first renormalon pole are determined. It 
is a measure for the ambiguity of the perturbative calculation and is 



assumed to dominate the power corrections. The corrections Ag% and 
c~| , A$5 are given as functions of the Bjorken variable x and turn out to be 

negligable. The anomalous dimensions of the leading twist operators 
are obtained in the next-to-leading order. 



It is well known that the perturbation series for moments of twist- 2 struc- 
ture functions is an asymptotic one. This property can be studied in detail 
in the 1/A^-limit, in which the complete series can be calculated explicitely. 
Formally this series is given by an integral over the positive real axis in the 
Borel plane. This integral is ambiguous due to singularities on the integra- 
tion path, the so called IR-renormalon poles. The residues of these poles are 
a measure for the ambiguity of the perturbative series. The so-called hypoth- 
esis of UV-dominance allows furthermore to interpret this ambiguity as an 
estimate for the power corrections, the program just sketched was already 
applied to all twist-2 structure functions except g^(x, Q 2 ) and g^x.Q 2 ) |3J 
- H. In this contribution we investigate these remaining two cases. Good 
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experimental data for power corrections to structure functions exists so far 
only for Fz,. The renormalon prediction fits this data surprisingly well. Let 
us note that similar renormalon analyses have recently been applied to a 
large range of other QCD observables ||] - [fTTfl . 

The Borel transformation of a perturbative series 

oo 

R = r^a + rid 2 + r2<2 3 = . . . = V] r n a n+1 , a = a s ■ Aix (1) 

n=0 

is defined as 

= (2) 

n=0 

R can be reobtained from its Borel transform by an integration over the 
positive real axis as 

r°° i 

R= / due~ u/a B [R}(u) . (3) 



o 



The coefficients of the original power series can also be obtained individually 
by taking the derivatives with respect to u 



k 



(4) 

u=0 



B[R] has pole singularities on the real axis, the so-called renormalons [12 
The poles on the positive w-axis, which are called IR-renormalons because 
they can be traced back to low momentum contribution to the loop integrals, 
lead to ambiguities in the retransformation @ because it is unclear wether 
they have to be passed above or below. The fact that no unambiguous 
retransformation exists reflects the fact that the perturbative expansions are 
asymptotic O, [Ljl . The ambiguities are of the order of magnitude 



AR = e- u/a Res(B[R}{ 



u 



(5) 

u=pole position 



and can be interpreted as a measure for generic uncertainties of perturbative 
predictions or in other words as an estimation for corrections beyond leading 
twist perturbation theory [4j] . 

In connection with the investigation of renormalons the NNA-approxi- 
mation (naive non-Abelianization) |15j is of particular interest because in the 
Borel it leads plane to an effective gluon propagator of a very simple form 
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allowing a calculation to all orders in the coupling constant. In the NNA- 
approximation we start with a restriction to the leading 1/iVy-terms (Nf-. 
number of quark flavors), which is the sum of all diagrams with only one 
exchanged gluon but an arbitrary number of quark loops. The missing terms 
are then approximated by the replacement Nf — 24it 2 (3q = Nf — 33/2, 
which corresponds to a restriction to the leading terms of an expansion in 
the one loop /3-function of QCD (3q = ^p-(ll — §-/V/). The resummation 
of all corresponding diagrams leads to the Borel transformed effective gluon 
propagator 

B[ 9 2 D^(k))(s) = (6) 
with the new variable s := f3$u [|5], ||, plj. c is a renormalization scheme 



dependent constant, in the MS-scheme c = —5/3. The expression (||) differs 
from the original gluon propagator essentially only by the power of k 2 in the 
denominator. Consequently a calculation of a Borel transform in the NNA- 
approximation in all orders of the coupling constant is not more complicated 
than the corresponding normal next-to-leading-order calculation. 

We now apply the described method to the structure functions g% and 
#5 measurable in polarized deep inelastic lepton-nucleon scattering. These 
structure functions are defined by the following terms in the decomposition 
of the hadronic scattering tensor 

m N (p li S u + SpPv) 2m N S-qg flu 

WV = 9s H 95 + ■■■ ■ (7) 

p ■ q p ■ q 



We adopted the conventions of [|1y[] , a comparison with other definitions used 
in the literature is given in ||18| . Since the contributions to W^ u shown in 
eq. (0) are parity violating they involve weak interactions. We are looking 
at the case of pure Z-boson exchange and the interference part of Z- and 7- 
exchange. In order to avoid operator mixing we consider the nonsinglet part, 
which is obtained by taking the difference between proton- and neutron- 
structure functions |1| [20| . To simplify the notation we write gj := g? — 
gj , j = 3,5. Neglecting higher twist contributions, the moments of the 
structure functions have the form 

9j, n ■= I dxx n g j (x,Q 2 ) = Aj, n C^ n {Q 2 ) , (8) 

J 

where Aj >n are the matrix elements of the leading twist nonsinglet operators 
and Cj tn (Q 2 ) the corresponding Wilson coefficients. The Wilson coefficients 
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can be calculated using their connection with the forward Compton scattering 
amplitude 

V = - m(P ^ + 2Ea 3 ,„C 3 ,n(Q 2 K +1 

P-Q n 

+ 2mg - g fr» 2 £ a 5 , n C 5 , n (g 2 )^ +1 + . . . , (9) 

P-Q n 

where t^ v and aj <n refer to quark states instead of nucleon states. Adopting 
a normalization where the non-vanishing Wilson coefficients take the form 

Cj, n (Q 2 )=l + 0(g 2 ) (10) 

the matrix elements of the leading twist operators are 

o 3 ,„ = 2 VA , n = 0,2,4... 

a 5 ,„ = Wl, n = 1,3,5... (11) 

with the vector coupling constant V and the axial coupling constant A. The 
Borel transformed Wilson coefficients are now obtained by the calculation of 
B[t^ u ] and comparing the result expanded in uj with eq. @. In the calcu- 
lations we have to handle the matrix 75 in d 7^ 4 dimensions. We use the 
t'Hooft-Veltman scheme 75 = i7°7 1 7 2 7 3 , {75,7^} = for \i = 0,1,2,3 and 
[75 j 7* 1 ] — otherwise We get 



d 



B[C 3 , n ](s) = C F £- 



xUd-2) 



q 2 J v Q 2 J (47r) d / 2 r(s + i)r(rf-i-s) 

T(s + n + 3 - d ^ 



2' 



nl 



+ |s + 2 _%,_4) ^r(. + » + 2)r(. + n + 4-S 



2/V v y / n!r(s + n + 4) 

d. nT(s + n + 2)T(s + n + 3-l) 

+-(6-d)- 



2 V ' n\T(s + n + A) 

+^(4(d - 4) 2 - 4(d - 4)(d - 2) + (d - 2) ; 
r(s + n + 2)r(s + n + 3- f) 

X 



n!r(s + n + 4) 
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+ (d-2-.)(rf-4)" r <' + ''+ 1)r( ' + " + 3 -' ) 
v A 7 n!_(s + n + 3) 

+4 ( s + 2 -^E r(s + " + 3 ~ f) 



27^0 A;!( S + A; + 1) 
_(s + A; + 3-f) 



' A;!(s + A; + 2) 



+(( 2a +4-d)^i-d)5: 

A;_(s + A; + 2-f 



d - 4 ^ ^ fc r(g + fc + 2- 1) 

k\(s + k + l) 



k=0 

" fcr( s + fc + 2 - 1 ) 

+ ^ A;!( S + A; + 2) 

+ (^ - 2 - *) ( 2 ( rf - 4 ) 2 - 2 ( rf - 4)(d - 2) + (d - 2) 2 ) 
" r(s + A; + 2-f) 

X t *!(* + * + !) 

x d-4 " A;r(s + A; + 2- |)1 , x 
+(d -4-2 S )( d -2- s) — g_L___.il (12) 



for n = 0, 2, 4 . . . and 



x{ (d-2) 



Q 2 J \ Q 2 J (4vr) d / 2 r(s+ l)r(d- 1-s) 
_(s + n + 3-f) 



n! 

d\/ ^r(s + n + l)r(s + n + 4- f) 



n!_(s + n + 3) 



+ ( s + 2--)(2-(d-4)): 

+ (, + 2 - g) ( ( d - 4) - 2) nV ^ s + n + ; 1)r(g + " + 3 " i} 

V 2/V v ; ; n ir( s + n + 3) 

+ ^2(d-2-s)(d - 4) 2 + -(6s+12-5d)(d - 4)(d - 2) 

+i(3d-8-4s)(d-2) 2 ) 



x £(s + n + l)r(s + n + 3- f) 



n\T(s + n + 3) 
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+ 4 a +2-S £ 



2/ fc f *:!(* + fc + 1) 
" kT(s + k + 2-l) 
t Q k\(s + k + l) 

^ - 2 - s) (2(d - 4) 2 - 2(d - 4)(d - 2) + (d - 2) ; 
xi:^^ 2 "^) (13) 



fc=0 



A;!(s + k + 1) 



for n = 1, 3, 5 . . . . 

Since the NNA-approximation is exact in one loop order we get the next- 
to-leading-order result from eqs. ([12]) and (|13D by taking s = according to 
eq. (U). An expansion in e = 2 — | leads to 

/ 4 4 4 

x - 4 + + + AS r 

V n + 1 n + 2 n + 3 

3 9 6 / 2 4 



2 + n + l n+3 + V + n+2 n + 3/ ™ 

+4 £^ +2 £ (*+4+2) M' (14 » 

C 5 ,„ = 1 + C ^{(I- 7 + ln ^)(-3 + 4- + 4- + 4 Sr 
(47r) 2 [Ve fi 2 / \ n + 1 n + 2 

4 6/2 2 \ 

"I + — (3 — + )S n 



n + 1 n+2 V n+1 n+2 



fc: 



-fc+1 



where S n is defined by if> n := J2k=i f • From the last two equations we read 
of the renormalization constants for the corresponding composite operators 
(defined by O r = Z^Oq) in the MS-scheme @. 

Z 95>n = 1 + ^^(7-7 + ln47r) 

/ 4 4 4 \ 

x(-4 + — T + — <+45„ , 16 
V n + 1 n + 2 ri + 3 / 
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1 + C r 



(4tt) 



7 + In Air J 
x ( - 3 + — r + — r + 4 S r 



n+1 n+2 



Finally we get for the anomalous dimensions 7 
one loop order 



Z dfj, 



9o 



(see e. g. 



(17) 



in 



7g3,n 



(4vr) 2 
.2 



_^ ^ + ^ 8T-1 

n + 1 n + 2 n + 3 ^ ' 



.9" 



(4vr)- 



6- 



4 



4 n 1 



n + 1 n + 2 A; 



(18) 
(19) 



To our knowledge these anomalous dimensions have not been calculated be- 
fore. Higher order results could be obtained in the — no longer exact — NNA- 
approximation as well using eq. (H). 

To investigate the renormalons we can set d = 4. From eqs. ( |I~2"D and ( |I~3"| ) 

we get 



C f 



' 1 — r\ s 

lie \ 



Q 2 ) (47r) 2 r(s + l)(2-s)s ' s-1 



J 2 £( i± n ± l)_ 2s 



r(s + n + 2) 



n\T(s + n + 4) 
^nT(s + n + 2)r(a + n + 1) ^(s + n + 2)r(s + n + 1) 

ra!r(s + n + 4) 



n\Y(s + n + 4) 



+2E 



^ [ 2S fc!(s + A; + 1) \l(s + fc + 2) Z k\(s + k + l) 

(20) 



r( g + fc + 1) _ r(s + fc + i) _ fcr(g + fc) 



o fcr( 8 + fc) 2g r(s + fc) 



'fc!(s + fc + 2) £;!(s + £; + !) 



£[C 5 ,n]0) = CV 



Q 2 / (47r) 2 r(s + l)(2-s)s 's-1 



x{2 



r(g + n + l) | ^ T(s + n + l) 

n\ n!(s + n + 2) 

2 



-2s- 



i(T(s + n + 1) 
nT(s + n + 3) 



+ 4(1 



(r(s + n + l) 
n!r(s + n + 3) 



+4E 

k=0 



' T(s + k + l) kT(s + k) T(s + k) 

S-TT7 ; 2—7 : — — S- 



k\(s + k + l) k\(s + k + l) k\(s + k + l) 



,(21) 



The pole at s = corresponds to the usual 1/e pole in dimensional regular- 
ization. In both cases we find two IR-renormalons for s = 1 and s = 2. The 
corrections corresponding to these renormalons are suppressed by factors 
or (^2) 2 respectively, which leads to the hypothesis that they should domi- 
nate these power corrections [fj, |3|, |J. For the dominant pole at s — 1 and 
taking fi 2 = Q 2 we find for the residues 



Res(B[C 3 , n ](s] 



Res(B[C 5tn ]( S ] 



C F e~ c 



2n + 6 



n 



n + 2 



16 24 
+ 



n + 3 n + 4 



(22) 



(47T) 



2n + 10 - 



n + 1 



n ] _ 



n+2 n+3 



(23) 



which are connected with the renormalon contributions according to eq. ([|) 
by 



AC^Q 2 



± 



'A 2 



Q 2 fa 



^Res( J B[Q, n ]( S 



+ 



(24) 



The sign of the corrections remains unknown. Correspondingly we get ac- 
cording to eq. (H) for the complete structure functions 



N 

gUQ 2 ) = A hn [T,cUQ 2 )(9 2 ) k + AC Jin (Q 2 



k=0 



(25) 



2\k 



with the perturbative expansion of the Wilson coefficients Cj tn = J2kCj t1l (g 
and for the renormalon corrections of the same structure functions 

Ag^ n (Q 2 ) = A hn /\C hn {Q 2 ) . (26) 

The unknown matrix elements A^ n are eliminated taking the ratio 



AC^Q 2 



Ekl CU9 2 ) k + A C hn (Q 2 
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± 



(f) iM B ^ 8 ))Li +0 (^) 



1 + Ek=i c jAQ 2 ) h (9 2 ) k + AC iin (Q 2 



^1 -1 



± 



'A 2 ' 



lRes(s[C7 i|n ]( a ; 



So in leading order the corrections are given as 

/A 2 \ 1 



+ 



s=l 



s=l 



The deternination of all moments is equivalent to expressing Agj\ 
convolution 

where the functions Cj R1 (y) defined by 

Res(B[C j , n ](s)) =~ 1 

s=l 



dyy n C m (y) 



(27) 



(28) 



X ) clS cL 



(29) 



(30) 



are obtained from eqs. ( p0[) and (|2T|): 



Ci m (y) 
Cl R \y)- 



(4vr) 2 
C F e~ c 



25'{y-l)-65{y-l)+A+Ay+16y 2 -2Ay 3 



25\y-l)-W5(y-l) + 8+4y+8y 2 



(i-y) 



(4tt) 2 { " > - " <" * (l-yJ+J' 

wner e is defined by /J dy /(y) ^q^ = /p 1 dy /(y ^ (1) ■ We 

quark distributions given in j24j and the parton model expressions 

2^<#<?i(Ag-Ag) , 



2^ 

#5 



ft* 



2x]Te^(Ag-Ag) 



-|,(3D 
(32) 
use the 

(33) 

(34) 
(35) 



2xgf 



95 



(36) 
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Figure 1: The fit for (full line) 
and the corresponding renor- 
malon contribution multiplied by 
a factor 10 (dotted line). The 
dashed lines show the size of the 
ambiguity for gf , i. e. #f ± Agf . 



Figure 2: The same as fig. [T] for 
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Figure 3: The same as fig. [T] for Figure 4: The same as fig. [3] for 



9s ■ 



95 ■ 
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We choose the momentum transfer to be Q 2 = 4GeV 2 . The integrals in eq. 
fl29p are evaluated numerically and the results are plotted in the figures [1] to 

We have thus completed our analysis of the renormalon ambiguities for 
all twist-2 structure functions. 
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